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Abstract 



Much has been studied on the spreading speed and traveling wave solutions for co- 
operative reaction-diffusion systems. In this paper, we shall establish the spreading 
speed for a large class of non-cooperative reaction-diffusion systems and character- 
ize the spreading speed as the slowest speed of a family of non-constant traveling 
wave solutions. As an application, our results are applied to a partially cooperative 
system describing interactions between ungulates and grass. 
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1 Introduction 



Fisher [TT] studied the nonlinear parabolic equation 

w t = w xx + w(l-w). (1.1) 

for the spatial spread of an advantageous gene in a population and conjec- 
tured c* is the asymptotic speed of propagation of the advantageous gene. His 
results show that (11. ip has a traveling wave solution of the form w(x + ct) if 
only if \c\ > c* = 2. Kolmogorov, Petrowski, and Piscounov [IE] proved the 
similar results with more general model. Those pioneering work along with 
the paper by Aronson and Weinberger [i~|2] confirmed the conjecture of Fisher 
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and established the speeding spreads for nonlinear parabolic equations. Lui 
[20] established the theory of spreading speeds for cooperative recursion sys- 
tems. In a series of papers, Weinberger, Lewis and Li [T7jl8f32|l3"3"] studied 
spreading speeds and traveling waves for more general cooperative recursion 
systems, and in particular, for quite general cooperative reaction-diffusion sys- 
tems by analyzing of traveling waves and the convergence of initial data to 
wave solutions. 

However, due to various biological or physical constrains, many reaction- 
diffusion systems are not necessarily cooperative. Thieme [27] showed that 
asymptotic spreading speed of a model with nonmonotone growth functions 
can still be obtained by constructing monotone functions. Weinberger, Kawasaki 
and Shigesada [36] discussed the reaction-diffusion model 

du\ , . r . , 

— — = di/Xui + ui[—a - ou\ + riu 2 \ 

du 2 (L2) 

— j- = d 2 Au 2 + u 2 r 2 [l -u 2 + h{ui)] 

where d%, a, S, r 1; d 2 ,r 2 are all positive parameters. This system describes the 
interaction between ungulates with linear density U\(x, t) and grass with linear 
density u 2 (x,t). The function h(ui) models the increase in the specific growth 
rate of the grass due to the presence of ungulates. When the density U\ is small 
the net effect of ungulates is increasingly beneficial, but as the density increases 
above a certain value, the benefits decrease with increasing U\. (11.21) is partially 
cooperative 2-species reaction-diffusion model, meaning that it is cooperative 
for small population densities but not for large ones. By employing comparison 
methods [36] established spreading speeds for (11.21) . In Section [5], we take the 
non- monotone Ricker function u\e~ Ul as h(u\), which is simpler than that of 
[36] , and apply our main theorem (Theorem 12.11) to (II. 2p . The application 
of our general theorem allows us to characterize the spreading speed as the 
slowest speed of traveling wave solutions to (ll.2p . which is new and was not 
proved in [36] . Non-cooperative react ion- diffusion systems frequently occur in 
other biological systems such as epidermal wound healing (see, Sherratt and 
Murray [25.26], Dale, Maini, Sherratt (7J). In a recent paper by the author [31] . 
spreading speeds and traveling waves for a non-cooperative reaction-diffusion 
model of epidermal wound healing were established. 



For related non-monotone integro-difference equations, Hsu and Zhao [M], Li, 
Lewis and Weinberger [19] extended the theory of spreading speed and estab- 
lished the existence of travel wave solutions. The author and Castillo-Chavez 
[3D] prove that a large class of nonmonotone integro-difference systems have 
spreading speeds and traveling wave solutions. Such an extension is largely 
based on the construction of two monotone operators with appropriate prop- 
erties and fixed point theorems in Banach spaces. A similar method was also 
used in Ma [23] and the author [2S] to prove the existence of traveling wave 



solutions of nonmonotone reaction-diffusion equations. 

In this paper, we shall establish the spreading speed for a general non-cooperative 
system (II .3p and characterize its spread speed as the slowest speed of a family 
of non-constant traveling wave solutions of ( II. 3p . Our main theorem (Theorem 
EH]) will be applied to ([L2D in Section ED 

We begin with some notation. We shall use i?, fc, fc 1 * 1 , /, f ± ) r ) u ) v to denote 
vectors in M, N or iV-vector valued functions , and x, y, £ the single variable in 
R. Let u = (Ui),v = (vi) G M. N , we write u > v if Ui > Vi for all i; and u ^> v if 

Ui > Vi for all i. A vector u is positive if u 3> 0. For any r = {ji) >0,rG R^ 
let 

[0,r] = {u:0<M<r,!iG R N } C M^ 

and 

C r = { u = (ui) : ^ e C(R, R), < Wi(x) < r i; x G R, z = 1, ..., A^}, 

where C(R, R) is the set of all continuous functions from R to R. 

Consider the system of react ion- diffusion equations 

u t = Du xx + f(u),xeR, t>0. (1.3) 

with 

u(x,0) =ito(aO, ^ R . (1-4) 

where w = (ui), D = diag(di, d,2, •••, dp?), di > 0, i = 1, ..., iV 

/( M ) = (/i( M )> /2(w), •••, /v(«)), 

i*o (x) is a bounded uniformly continuous function on R. In this paper, by 
a solution we mean a twice continuously differentiable function u satisfying 
appropriate equation in R and an initial condition. 

In order to deal with non-cooperative system, we shall assume that there are 
additional two monotone operators f ± , one lies above and another below / 
with the corresponding equations 

ut = Du xx + f + (u), x G R, t > 0. (1.5) 

Ut = Du xx + f-{u),xeR,t>0. (1.6) 

Such an assumption will enable us to make use of the corresponding results 
for cooperative systems in [20,32] to establish spreading speeds for (11.31) . 

(HI) i. Let k + = (£;+) >> and / : [0, k + ] ->■ R N be a continuous and 
twice piecewise continuous differentiable function. Assume that C&+ is 
an invariant set of (II. 3p in the sense that for any given u G C&+, the 



solution of (II. 3p with the initial condition u® exists and remains in Ck+ 
for t G [0,oo). 
ii. Let << k~ = (&j~) < k = (ki) < k + . Assume there exist continu- 
ous and twice piecewise continuous differentiable function f ± = (f^) : 
[0, k + ] -> R N such that for u G [0, fc+] 

/-(«) < /(«) < /+(«). 

iii. /(0) = /(A;) = and there is no other positive equilibrium of / between 
and k. f ± (0) = f ± (k ± ) = 0. There is no other positive equilibrium of 
f^ 1 between and k . 

iv. (|1.5p and (11.61) are cooperative (i.e. diff(u) > for u G [0, fc ± ],i 7^ j). 
f ± (u), f(u) have the same Jacobian matrix /'(0) at u — 0. 

A traveling wave solution u of (II .3p is a solution of the form u = u(x + ct) , u G 
1 N ) . Substituting u(x, t) = u(x + ct) into (11.31) and letting £ = x + ct, 



we obtain the wave equation 

Du"(0 - «/(0 + /( u (0) = 0, (el. (1.7) 

Now if we look for a solution of the form (uj) = (e x ^r] l x J , A > 0, ^a = (ry A ) >> 
for the linearization of (11.71) at the origin, we arrive at the following system 
equation 

diag(^A 2 - c\) Vx + f'(0)rj x = 

which can be rewritten as the following eigenvalue problem 

1 , . , 

-^A^A = c?7a, (1.8) 

where 

A x = (a A J ) = diag(^A 2 ) + /'(0) 

The matrix /'(0) has nonnegative off diagonal elements. In fact, there is a 
constant a such that f'(0)+al has nonnegative entries, where I is the identity 
matrix. 

By reordering the coordinates, we can assume that /'(0) is in block lower tri- 
angular form, in which all the diagonal blocks are irreducible or 1 by 1 zero 
matrix. A matrix is irreducible if it is not similar to a lower triangular block 
matrix with two blocks via a permutation. From the Perron- Frobenius theo- 
rem any irreducible matrix A with nonnegative entries has a unique principal 
positive eigenvalue( which is the spectral radius of A, p(A)) with a correspond- 
ing principal eigenvector with strictly positive coordinates. For an irreducible 
matrix A with nonnegative off diagonal elements, we shall call the eigenvalue 
p(A + al) — a of A, which has the same positive eigenvector, the principal 




Fig. 1. The red curve is ^(A). The minimum of $(A) is c* . For c > c*, the left 
solution of ^(A) = c is A c . 

eigenvalue of A (see e.g. [I3B2])- Let 

■qr(A) = p(A + aI) -a. 

Here A + al is irreducible and nonnegative, and p(A + al) is the spectral 
radius of A + al. 

We shall need the following assumption (H2). Notice that (H2) is assumed for 
A = in [32]. However, with (H2), we are able to obtain better estimates for 
traveling solutions and the minimum speed c*, see Lemma [TTT1 As a result, for 
the example in Section [51 a slightly stronger assumption (d\ > g^) than [36] 
is be assumed. 



(H2) Assume that, for each A > 0, A\ is in block lower triangular form and 
the first diagonal block has the positive principal eigenvalue ty(A\), 
and ty(A\) is strictly larger than the principal eigenvalues of all other 
diagonal blocks, and that there is a positive eigenvector v\ = ( 



of Ax corresponding to ^(Ax) 
with respect to A. 



*4) » o 

And further assume that vx is continuous 



Let 



$(A) = j*(A x ) 



>0. 



According to Lemma 11.11 we can expect the graph of $ as in Fig. [TJ For the 
example in Section \5\ $ is a strictly convex function of A and, clearly satisfies 
Lemma 11.11 



Now we state Lemma fl.lt which shall enable us to give accurate asymptotic 
estimates of traveling solutions. Lemma [TTT1 is a analogous result in Weinberger 
[35] and Lui [20]. However, due to the fact that f'(0) is only quasi-positive and 
the elements of Ax are not necessarily log convex, some of its proofs here are 
different from Lui [20] . A theorem on the convexity of the dominant eigenvalue 
of matrices due to Cohen [5] is used to show that ^(Ax) is convex function of 



A. 

There are two direct consequences of Lemma 11.11 First one, it improves 
Theorem 4.2], which will be used in this paper, by eliminating the case (b) in 
[32| Theorem 4.2] because of Lemma ll.ll (7). Second, Lemma lTTTl (8) will allow 
us to construct explicit lower solutions and therefore asymptotic behavior of 
traveling solutions of (II. 3ft can be obtained. 

Lemma 1.1 Assume that {HI) — (H2) hold. Then 

(1) $(A) -> oo as A -> 0; 

(2) $(A) -¥ oo as A -» oo; 

(%) $(A) zs decreasing as A = + ; 

(^ \I/(A\) zs a convex function of X > 0; 

ftfj 3>'(A) changes sign at most once on (0, oo) 

(^7j $(A) /ias £/ie minimum 

c* = inf $(A) > 

A>0 V y 

at a finite X. 
(8) For each c > c* , there exist A c > and 7 e (1, 2) such that 



That is 



and 



$(A C ) = c, $(7A C ) < c. 

t-^4a c ^a c = $(A c )z/ Ac = cv Ac 
-A 7Ac z/ 7Ac = $(7A c )z/ 7Ac < cz/ 7Ac 



7 A C 

where u Ac » 0, z/ 7 a c >> are positive eigenvectors of j-A/ ili ,-^-A 7 / ic 
corresponding to eigenvalues 3>(A C ) ana 7 $(7A C ) respectively. 



PROOF. We only need to prove those different from [35.J20J. The proof of the 
convexity of ^>(A\) is similar to that in Crooks [6] for matrices with positive 
off-diagonal elements. It is easily seen that ty(A\) = p(A\ + al) — a is non- 
decreasing function of A > ([131 Theorem 8.1.18]). Further, a theorem on 
the convexity of the dominant eigenvalue of matrices due to Cohen [5] states 
that for any positive diagonal matrices D 1; D 2 and t G (0, 1), 

q,(t Dl + (1 + t )D 2 + /'(0)) < t^{D 1 + /'(0)) + (1 - t)¥(D 2 + /'(0)) 

As before ^(A) is the principle eigenvalue of A. Now if a\,a 2 £ K and t G 
(0,1), 

(toi + (1 - t)a 2 f < ta\ + (1 - t)«2- 



This implies that 

*(A tAl+(1 _ t)A2 ) = *((tAi + (1 - t)A 2 ) 2 D + /'(O)) 
<*(£\?£>+(l-t)A§£> + /'(())) 

+ a - tm%D + no)) 

= t*(A Al ) + (l-t)*(A Aa ) 

Since \I/(A A ) is a simple root of the characteristic equation of an irreducible 
block, it can be shown that ty(A\) is twice continuously differentiable on R. 
Thus 

*"(A) > 



and a calculation shows 



[A$(A)]' = #'(A) 
$'(A) = -^'(A)-$(A)] 



and 



(A 2 $'(A))' = A*"(A) > 0. 

(6) is a consequence of the above inequalities. As for (2), we need to prove 
that lirn A _ ) . 0O ^ a a) = oo. In fact, there exists an e > such that all diagonal 
elements of D — el are strictly positive, then ty(D — el) > and choose A 
large enough so that 

= *((Z?- e I)A 2 + (eA 2 / + /'(()))) 

>ty((D-eI)\ 2 ) 
= \ 2 ^(D-eI^ 



Thus lim A ^oo A = oo. As we discussed before, (H2) implies the existence of 
positive eigenvector v\ » corresponding to $(A A ). The first statement of 
(8) is a consequence of (l)-(7). The second statement of (8) is just a rephrase 
of the fact that v\ » is a eigenvector of jA\ corresponding to eigenvalue 
$(4\) for A = A c and 7A C . 

D 



In addition to (H1-H2), we also need assumption (H3) which only requires 
the nonlinearity grow less than its linearization along the particular function 
v\t~ Xx [32]. Such a condition can be satisfied for many biological systems. 

(H3) Assume that for any a > 0, A > 

/ ± (a^A) < af'(0)u x , where v x = (u{). 



We now recall results on the spreading speeds in Weinberger, Lewis and Li 
[32] and Lui [20]. While Theorem 4.1 [32] holds for non cooperative reaction- 
diffusion systems, it does require that the reaction-diffusion system has a sin- 
gle speed. In general, such a condition is very difficult to verify. In the same 
section, for cooperative systems, Theorem 4.2 in [32J provides sufficient con- 
ditions to have a single speed. The following theorem combines the results of 
Theorems 4.1 and 4.2 in [32], which can be a consequence of Theorems 3.1 
and 3.2 for discrete-time recursions in Lui |20|. 



Theorem 1.2 (Weinberger, Lewis and Li {32}) Assume (HI) — (H3) hold and 
U.3)) is cooperative. Then the following statements are valid: 

(i) For any u G Ck with compact support, let u(x,t) be the solution of 111.3)) 
with fl-Mp . Then 

lim sup u(x, t) = 0, for c > c* 

t ~^° c \x\>ct 

(ii) For any strictly positive vector u> G M. N , there is a positive R^ with the 
property that if Uq G C& and Uq > u on an interval of length 2R UJ , then 
the solution u[x,t) of M.3)) with U.4\ ) satisfies 

liminf inf u(x,t) = k, for < c < c* 

«->oo \ x \<tc 

In another paper [18] . for cooperative systems, Li, Weinberger and Lewis 
established that the slowest spreading speed c* can always be characterized 
as the slowest speed of a family of traveling waves. These results describe 
the properties of spreading speed c* for monotone systems. Based on these 
spreading results for cooperative systems, we will discuss analogous spreading 
speed results for non cooperative systems. 



2 Main Results 



Our new contributions in this paper are to establish the spreading speed (The- 
orem [5]T] (i-ii)) for general non-cooperative reaction-diffusion systems ( II. 3p . 
and further characterize the spreading speed as the speed of the slowest non- 
constant traveling wave solutions (Theorem 12.11 (iii-v)). 

Although the existence of traveling wave solutions for cooperative systems 
are known (see,e.g. [IB]), we shall prove the existence of traveling wave so- 
lutions for both cooperative and non cooperative systems as our proofs for 
non cooperative systems are based on those for cooperative systems. Further, 
in additions to the existence of traveling wave solutions, we shall be able to 



obtain asymptotic behavior of the traveling wave solutions in terms of eigen- 
values and eigenvectors for both cooperative and non cooperative systems. 
The following theorem is our main results. 

Theorem 2.1 Assume (HI) — (H3) hold. Then the following statements are 
valid: 

(i.) For any u G Ck with compact support, the solution u(x,t) of M.3\) with 
\l-4\) satisfies 

lim sup u(x, t) = 0, for c > c* 

*-*-°° \x\>tc 

(ii.) For any vector uj G M. n ,uj >> 0, there is a positive R w with the property 
that ifu G Ck and u > uj on an interval of length 1R^, then the solution 
u(x, t) of M.3\) with \1.4\) satisfies 

k~ < liminf inf u(x,t) < k + , for < c < c* 

- t^oo \ x \< tc V ' - 

(Hi.) For each c > c* U.ty) admits a traveling wave solution u = u(x + ct) such 
thatO «m(0 < fc + ,(el, 

k~ < liminf it(£) < limsupu(£) < k + 

and 

lim w(0e" Ac€ = v Kc . (2.9) 

£->-oo 

If, in addition, U.fy) is cooperative in Ck, then u is nondecreasing on R. 
(iv.) For c = c* U.3\) admits a nonconstant traveling wave solution u = u(x + 
ct) such that < u{£) <k + ,£e R, 

kT < hminf u(£) < limsupw(£) < k + . 

(v.) ForO < c < c* M.3\) does not admit a traveling wave solution u = u(x+ct) 
with liminf^_ >00 'u(^) >> and u(— 00) = 0. 

Remark 2.2 In many cases, f ± can be taken as piecewise functions consisting 
of / and appropriate constants as demonstrated in Section In order to have 
a better estimate for the traveling wave solution u for non cooperative systems, 
it is desirable to choose two function f ± which are close enough. The smallest 
monotone function above / and the largest monotone function below / are 
natural choices of f^ if they satisfy other requirements, See [2T|14|I19] for the 
discussion for scalar cases and [SB] for a partially cooperative reaction-diffusion 
system. Our construct of /" in Section is different from the previous papers. 

Remark 2.3 The invariant set of (HI) (i) can often be established by Com- 
parison Principle 13.11 In fact, for a given u G Ck+, let u(x,t) be the solution 



of (11. 3p with the initial condition uq. If we can choose appropriate / , f + so 
that f-(u) < f{u) < f + (u) for all u e R N , it follows that 

# " ^*4 " / + (^ + ) = = «t - Du** - /(«) > Ut - Du xx - f + {u). 

and 

- DO - /-(0) = = ^ - Du^ - /(«) <u t - Du xx - f-(u). 

Comparison Principle 13.11 implies that 

<u(x,t) <k + ,xeR,t>0. 

Now according to Smoller [2U Theorem 14.4] (I1.3P (and also (jl.5p . ( 11.61) ) has 
a solution u for £ e [0, cxd) and < w < /c + if the initial value u is uniformly 
continuous on R. Now we can establish an invariance set for (I5.26P in Section 
First we can extend h, h + in (I5.26P and (15.321) to zero for W\ < 0. Now let 
/ + be the reaction terms in (I5.32p . And further let / _ be the reaction terms 
in (I5.26P with h being replaced by the constant zero function for all w\ e R. 
From the above discussion, we can see that [0, k + ] is an invariance set for 

Remark 2.4 When (11.31) is cooperative in Ck, then f^ 1 = f. 

Remark 2.5 As indicated in [32J, if / is not defined everywhere, (H3) can 
be replaced by (H3'). Without extra assumptions, (H3') can be verified with 
a slightly complicated computation. 

H3' For each A > 0, let v ± = (minjfcj , au\\). Assume that for any a > 

/^(x)) < a/'(0)K). 

The paper is organized as follows. Theorem 12.11 (i)-(ii) shall be proved in 
Section [3] and Theorem 12.11 (iii)-(v) in Section HI 



3 The Spreading Speed 

3.1 Comparison Principle 



We state the following comparison theorem for cooperative systems of reaction- 
diffusion equations in Weinberger, Kawasaki and Shigesada [36J . The com- 
parison principle is a consequence of the maximum principle (see, e.g., Protter 
and Weingberger 



10 



Theorem 3.1 Let D be a positive definite diagonal matrix. Assume F = (Fj) 
is vector-valued functions in M. N are continuous and piecewise continuously 
differentiable in R and the underling system is cooperative in the sense that 
for each j , Fj is nondecreasing in all but the jth component. Suppose that 
u(x,t),v(x,t) satisfy 

ut - Du xx - F(u) <v t - Dv xx - F(v) (3.10) 

Ifu(x,t ) < v(x,t ), x G R, then 

u(x,t) < v(x,t), x G R, t > t . 

We are now able to prove Parts (i) and (ii) of Theorem 12.11 
3.2 Proof of Parts (i) and (ii) of Theorem \2.1\ 



Part (i). For a given uq G Ck with compact support, let u + (x, t) be the solutions 
of (II. 5p with the same initial condition w , then Comparison Principle 13.11 
implies that u + (x,t) G C&+ and 



u+ - Dut x - f + (u + ) = Q = ut- Du xx - f(u) >u t - Du xx - f + (u). 
and 

- DO - f~ (0) = = u t - Du xx - f(u) <u t - Du xx - /"(«). 
Comparison Principle 13.11 further implies that 

< u(x,t) < u + (x,t),x G R,t > 0. 
Thus for any c > c*, it follows from Theorem 11.21 (i) that 



Jim sup u + (x,t) = 0, 
and hence 



i_> °° \x\>tc 



lim sup u(x,t) = 0, 

^°°|x|>te 

Part (ii). According to Theorem 11.21 (ii). for any strictly positive constant u>, 
there is a positive R^ (choose the larger one between the R^ for (jl.5p and the 
Rui for (11. 6p ) with the property that if u > u on an interval of length 2R U1 , 
then the solutions m ± (x, t) of (II. 5p and (jl.6p with the same initial value u are 



11 



in Ck+ and satisfy 



liminf inf u (x,t) — k , for < c < c*. 

t^oo \ x \<tc 

As before we have 

u? - Du+ X - f + {u + ) = = u t - Du xx - f(u) >u t - Du xx - f + (u) 
and 

u~ - Du~ x - f~(u~) = = ut - Du xx - f(u) <u t - Du xx - f~{u). 

Thus, Comparison Principle 13.11 implies that 

u~(x,t) < u(x,t) < u + (x,t),x 6l,i>0. 
Thus for any c < c* , it follow from Theorem 11.21 (ii) that 

liminf inf u^ixA) = k , 

t-^oo \x\<ct 

and hence 

k~ < liminf inf u(x, t) < k + . 

— t^oo \ x \<ct ~ 

a 



4 The characterization of c* as the slowest speeds of traveling waves 

4-1 Equivalent integral equations and their upper and lower solutions 



In order to establish the existence of travel wave solutions, we fist set up 
equivalent integral equations. Similar equivalent integral equations were also 
used before, see for example, Wu and Zou [38], Ma [22|23] and the author [29] . 
For the convenience of analysis, in this paper and [29], both Aij, A2J are chosen 
to be positive, and —Xu, \2i are solutions of ( 14. lip . 

Let (3 > max{\difj(x)\,x e [0, k + ],i,j = 1,...,N] > 0. For c > c*, the two 
solutions of the following equations, 

di\ 2 -c\-[3 = 0,i = 1,...,N (4.11) 

are —Xu and A 2 j where 

Xu ~ u > °' A2J " u > °- 



12 



We choose (3 sufficiently large so that 

X 2i >X li >2A c ,i = l,...,N. (4.12) 

Let u = (ui) G Ck and define a operator 7~[u] = (7i[u\) by 



U>i{Mi "T Mi) \Joo 



(4.13) 



where 

H l (u(s))=Pu l (s) + f l (u(s)), 

7i[u],i = 1,...,N is defined on K if Hi(u),i = 1,2 is a bounded continuous 
function. In fact, the following identity holds 



di{Xii + A2i) ^ J -co 






roc * 

Xl '^- s) (3ds + / e X2 ^- s) (3ds) 



g /I | 1\_ g (4.14) 

^i(Aii + A2i) ^Aii A2r ^i(AijA2j) 
= 1. 

We shall show that a fixed point u of T or solution of the equation 

w(0 = T*[u](0 ^K, (4.15) 

is a traveling wave solution of (11.31) in Lemma 14.11 
Lemma 4.1 Assume (HI — H2) hold. If u E Ck is a fixed point ofT[u], 

i/ien u is a solution of (1.1). 



PROOF. Note that Hi(u(s)) are continuous functions on M. Thus T[w](£) is 
defined and differentiable on R. Direct calculations show 

(7fa](fl)' = 7a ^1 J - *u T e^^M*))^ 

"U^li i" ^2i) V J-oo 

+ X 2i J °e x ^- s) H i (u(s))ds^ 



13 



and 

CM(O)" = 7nnrTr( A ii f e-^-^OO)* 

/■oo 

+ A|i / e^-^Ms))^ 
-Ai i fl- i («(0)-A WJ H i («(0) N 



Noting that — Aij,A 2i are solutions of (14.111) . one can evaluate the following 
expression 

(7Mfl)"-c(7-[«]«))'-/37-[«]«) 

= /37-M(0-^(m(0)-/37IM(0 
= -^«(u(0)- 

Now if u(f) = T[tt](0, f e K, then « and is a solution of ([117]). □ 

We now define upper and lower solutions of (14.151) . <p + and <p~ , which are only 
continuous on R. Similar upper and lower solutions have been frequently used 
in the literatures. See Diekmann [8], Weinberger [35], Liu [20], Weinberger, 
Lewis and Li [32], Rass and Radcliffe [21) . Weng and Zhao [37] and more 
recently, Ma [23], Fang and Zhao [9] and Wang [2U][30"] . In particular, it is 
believed that the vector- valued lower solutions of the form in this paper first 
appeared in [37] for multi-type SIS epidemic models. In this paper, the upper 
and lower solutions here are defined for general react ion- diffusion systems and 
we calculate the associated integrals to verify the validity of the upper and 
lower solutions. 

Definition 4.2 A bounded continuous function u = (itj) G C(R, [0, oo) ) is 
an upper solution of (14.151) if 

7I[m](0 < Ui{£), for all £ G R,i = 1, ..., N; 

a bounded continuous function u = («j) G C(R, [0, oo)^) is a lower solution 
of (1215]) if 

7~M(0 > u^), for all £ el,i= 1,...,JV. 

Let c > c* and consider the positive eigenvalue A c and corresponding eigen- 
vector z/a c = {p\ ) in Lemma [1.11 and 7 > 1, q > 1. Define 

+ (O = (tf), 
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Fig. 2. For each i, the curve above is (fif and the below is (\> i . 

where 

<PJ = min{k ij ui c e A °t},i= 1,...,N, £ G K; 

and 

<T(0 = (ft"). 

0" = max{0, v\e^ - qu^e^}, i = 1, ..., N, £ E R. 

hi hi- 

lt is clear that if £ > -^, 0+(£) = &, and £ < ^S <#(0 = <e A ^ = 
1,...,N. 

Similarly, if £ > (1 _^ e > & (0 = °> and for £ < ( i-^a c » 

0r(O = ^l c e AcC - ^ Ac e 7AcC , z = 1, -, N. 
We choose q > 1 large enough that 

ln(g^) IhjJl 

7 TT< —^,i = l,...,N 

(1- 7 )A C A c 
and then 

We now state that + and 0~ are upper and lower solution of (I4.15P respec- 
tively and their proofs will be carried out in Appendix through careful analysis 
of the associated integrals. 

Lemma 4.3 Assume (HI) — (H3) hold and U.3\) is cooperative. For any c > 
c* , <f) + defined above is an upper solution of fl^.i5p . 



Lemma 4.4 Assume (HI) — (H3) hold and 7 satisfies &6.44 ) '■ For any c > c* , 



defined above is a lower solution of Ij4-15\ ) if q (which is independent of £) 
is sufficiently large. 
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4-2 Proof of Theorem \2.1\ (Hi) when U.3\) is cooperative 



In this section, we assume that (jl.3p is cooperative and prove Theorem 12.11 
(iii). In this case, f^ 1 = f. Many results in this section are standard and for the 
verification of continuity and compactness of the operator. See, for example, 
Ma [2"2"f2"3"] and Wang [2H]- Define the following Banach space 



E Q = { u = ( Ui ) : Ui G C(R),sup|tti(f)|e~'* < oo,z = l,...,iV} 

£eR 



equipped with weighted norm 



N 

^sup|w;(0|e~^, 



where C(R) is the set of all continuous functions on R and g is a positive 
constant such that g < A c . It follows that <p + G £ e and 0~ G £ g . Consider the 
following set 

A = {u = (m) : Ui G C(R) G ^,07(0 <Ui< <#(£),£ e R,z = 1, ...,N.} 

We shall show the following lemma. 

Lemma 4.5 Assume {HI) — (H3) hold and difj > 0, i ^ j on [0, k]. Then T 
defined in {4-13^ is monotone and therefore T(A) Q A. Furthermore, Tl[u\ is 
nondecreasing if u G A and all of Ui are nondecreasing. 



PROOF. Note that Hi(u(£)) and T[w](£) are bounded continuous functions 
on Rif u G A. Note f3 > max{\difi(u)\,u G [Q,k],i = 1,-,N} > 0, $&•(«) > 

0,u G [0,/c],i 7^ j. For any u = (ui),v = (itj) G A with tti(£) > fj(£)>£ G R„ 
we have, for £ G R 



fli(tt(0) - #(«(£)) 

= 0MO - vM)) + jf 1 ^MO + (1 - s)v(0)ds(u(0 - v(0) (4-16) 
>0. 
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If u G A and Ui are nondecreasing, we get, for i = 1, ..., N, £ G R and £1 > 0, 



/■oo 

+ / e A2 ' K+ ^ s) iJi(!i(s))rfs 



e -^- s )#.( w ( s ))d s 



-oo 
oo 



e ^*«-»)fl i (u( 5 ))ds) 

'* 1 , /* ^ 4 - 17 ) 

ttt -r r (/ e-^«-.)fl i ( Tt ( a + ei )) ( i s 

e- Ali «- s ^i(w(s))ds 



-oo 

oo 






It follows from d4T6|) that 7l[u](£ + &) - Ti[u](£) > for £ G R and ft > 0. 

□ 



Now we shall show that T[u] is continuous and maps a bounded set in A into 
a compact set. 

Lemma 4.6 Assume (HI) — (i/3) /ioW. T/ien T : A — > £ e is continuous with 
the weighted norm ||.|| e . 



PROOF. Let 

L = maxildifjiu^u G [0,fc],i = 1,...,N}. 

For any it = (ui), v = (vi) G A, we have, for £ G R 

< /?k(0 - t^flle-* + | jf 1 ^MO + (1 - *)«(£))&(«(£) - «(0)|e" rf 

<(/3 + L)||«-t;||, 

(4.18) 
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Thus, we obtain 



|7IM(0-7M0|e^ 
< 7T I x A f e~^-s) lHMs)) _ HMs))lds 

{Mi "T A2i) x J-oo 

+ / e A " ({ ->|fr i («)(s) - ffi(i.)(s)|ds)e-« 5 

< <4 l >H"7"'( r .-*»«->«-*. (4 ' i9) 

(Xu + g)(\2i - g) (Aii + A 2 i) 



and 



Jlle - mm,{(A li H-^)(A 2i -^)} 11 



e- 



Thus, T[w] is continuous. □ 



Lemma 4.7 Assume (HI) — (H3) hold. Then the set T(A) is relatively com- 
pact in E e . 



PROOF. Let Mi = max^g^gK i7j («(£)) < oo,i = 1, ...,N. Recall that 



di(X 



- . [/ e - x ^- s Us+ e x ^ t - s Us} = ^. 

\i + A2i) J-oo Jt (j 



If u G A, £ G R and 5 > ( without loss of generality), we have, i = 1, ..., N 

T i [u]ti + S)-T i [um 
= \ ( f H5 e-^ +s ^H t (u(s))ds 

/•oo 

+ / e X2i ^ +s - a) Hi(u(s))ds 

J£+5 

e- Xli ^- s) Hi(u(s))ds 

— oo 
oo . 

e ^K«-»)fl.( u ( s ))ds] 
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(ii(Aii + A 






(4.20) 



2iJ \ "'-oo 

+ / e A2 ' K - s) (e A2 ^^(u(s)) - fli(«(s)))ds 
+ f^' e- x ^ +5 -^Hi(u(s))ds 

C+ *e A «« + *-)fl i (M(s))d3 



and 



7J[»]« + (!) - 7J[«](0| < m(K{|e- A » s - 1|, |e A » s - 1|}^ 

+ f .„^, ,+faW. "< 



^j(Aii + A2i) rfj(Aii + A2i) 

Thus we establish that 

]im(%[u](£ + 6) -7J[m](0) = 0, uniformly for all u G A £ G R, i = 1,...,N. 

8— >0 

(4.21) 
Take any sequence (u n ) = (uf) G A and let f™ = (t>™) = T[u n ]. From Lemma 
14.51 and (14.211) . (t> n ) is uniformly bounded on R and uniformly equicontinuous. 
For I m = [—m,m], m G N, by Ascoli's theorem and the standard diagonal 
process, we can construct subsequences {u nm ) of (u n ) such that there is a 
function v = (v j), fj G C(— oo, oo), i = 1, ..., JV and (v n,n = T[u nm ]) uniformly 
converges to v on each I m for m G N. Now we need to show that v E A and 
||^™ m - ^|| e -^ as n m -> oo. By Lemma 1431 4 r (O < ^(0 < 0i~(£)^ = 
1, ..., iV for all £ G R, and therefore v G A. Note that 

lim (0+(£) - 0"(£))e-^ = 0,2 = 1, ...,iV. 

£— >±oo 

For any e > 0, we can find Kq > such that if |£| > Kq, then, for all m G N 
K m (0 - v % \e^ < (0+(£) - 0"(£))e-^ < e,z = 1, ...,iV. 
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On the other hand, on [—I m ,I m ], (v nm ) uniformly converges to v. Thus there 



exists a L > such that, for n m > L 



K m (C)-Vi\e-^ < e, £ G [-K ,K ],i = 1,...,N. 
Consequently, if n m > L, the following inequality is true for all £ G 

\Vi m (Z)-v i \e-'*<e,i = l,...,N. 



Thus \\v nm - v\\ e ->■ as n m ->• oo. □ 



Now we are in a position to prove Theorem 12.11 when (11.31) is cooperative. 
Define the following iteration 

u 1 = K 1 ) = T[</>+], u n+1 = «) = 7>"], n > 1. (4.22) 



From Lemmas 14. 3| 14.41 14. 5| u n is nondecreasing on R and 

0"(O < < +1 (0 < <(0 < 0+(£U e R, n > l,n = 1, ...,iV. 

By Lemma |4~T1 and monotonicity of (« n ), there is u G ^4 such that lim n _ >00 ||'U n — 
u|| e = 0. Lemma [4.61 implies that T[u] = u. Furthermore, u is nondecreasing. 
It is clear that lini£_ ) ._ 00 'Uj(£) — 0, i — 1,...,N. Assume that lim^oo Wj(£) = 
hi, i — 1, ..., N k'i > 0, i — 1, ..., iV because of it G A Applying the Dominated 
convergence theorem to (14.13J) . we get k\ = \{fik\ + /j(/ei, ..., k' n ) By (H2), 
k[ = ki. Finally, note that 

"l(^ - <? e7Ac? ) < «i(0 < <e A ^,£ G R. 
We immediately obtain 

lim Ui (0e- A °t = u\ c ,i= 1,...,N. (4.23) 

This completes the proof of Theorem 12.11 (iii) when f |1.3|) is cooperative. 
^.5 Proof of Theorem \2.1\ (iii) 



PROOF. Theorem 12.11 (iii) is proved when fl 1 . 3 p is cooperative in the last 
section. Now we need to prove it in the general case. In order to find traveling 
waves for ( II. 3p . we will apply the Schauder's fixed point theorem. 

Let u = (ui) G A and define two integral operators 

7*M = (T^M) 
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for / and / + 

= TTT \, j / € e- A ««- S )F±(«(5))d S H-r e ^«- S )^(«(5))d5] (424) 

and 

Hi t (tt( S ))=/3u i ( fl ) + )f(« (fl )). 

As in Section I4T2"} both T + and T~ are monotone. In view of Section H~2l and 
the fact that / _ is nondecreasing, there exists a nondecreasing fixed point 
u~ = (u^) of T~ such that T~[w~] = u~, lim^ 00 u~(^) = fc~,z = 1,...,N, 
and lim^_ i ._ oo M7(0 = 0,i = 1,...,N. Furthermore, lim^_ 00 M~(^)e~ Ac ^ = 
v\ ,i = 1, ..., N. According to Lemma [4. 3[ <p + (with fc being replaced by k^) 
is also a upper solution of T ± because the proof of Lemma 14.31 is still valid if 
/ is replaced by f ± . Let 

f + (0 = (#(£)), 

where 

if(0=mm{kt,u l Ac e A ^}, i = 1,...,N, £ G R; 

It follows that Ui{£,) < (f)f, £ e R, i = 1, ..., N. Now let 

B = {u : u = {m) e £ ,u-{$ < Ui {$ < 0+(O,e e (-00,00), i = 1, ..., JV}, 

(4-25) 
where £g is defined in Section 14.21 It is clear that B is a bounded nonempty 
closed convex subset in E e . Furthermore, we have, for any u = (u{) G B 



U: 



Tr[u-} < %~[u] < %[u\ < % + [u] < % + [^} < #,i = 1,...,N. 



Therefore, T '■ B — >■ B. Note that the proof of Lemmas 14.61 14.71 is valid if 
(jl.3p is not cooperative. In the same way as in Lemmas 14. 6(14.71 , we can show 
that T : B — > B is continuous and maps bounded sets into compact sets. 
Therefore, the Schauder Fixed Point Theorem shows that the operator T has 
a fixed point u in B, which is a traveling wave solution of fll.3p for c > c*. 
Since u~(£,) < Ui(£) < $f(£)>£ e (— 00,00), i = 1,...,N, it is easy to see that 
for i = 1, ...,iV, lini£_ ) ._ 00 «*-(£) = 0, lim£_»_ 00 'U i (£)e' 



-A c £ 



IV 



A c 5 



fc < liminf u(£) < limsupu(£) < k + 



and < «-(£) < ^(£) < fc+ £ G (-00,00). D 
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4-4 Proof of Theorem \2.1\ (iv) 



PROOF. We adopt the limiting approach in [Ij to prove Theorem 12.11 (iv). 
For each n G N, choose c n > c* such that lin^^oo c n = c*. According to 
Theorem 12.11 (iii), for each c n there is a traveling wave solution u n = (u™) of 
( II. 3p such that 

U n = T[u n ](£). 

and 

k~ < liminf <(0 < limsup<(0 < kf,i= 1, ...,N. 

As it has shown in ( I4.20p . (u n ) is equicontinuous and uniformly bounded on 
R, the Ascoli's theorem implies that there is vector valued continuous function 
u = {u 1 ) on K. and subsequence (u nm ) of (u n ) such that 

uniformly in £ on any compact interval of R. Further in view of the dominated 
convergence theorem we have 

« = T[«](0. 

Here the underlying Ah, \2i of T is dependent on c and continuous functions 
of c. Thus u is a traveling solution of ( jl .31) for c = c*. Since, for each c n , u n G B 
where B is defined in (14.251) . it is easy to see that u satisfies 

kr < liminf w(£) < limsupu(£) < k + 

Because of the translation invariance of u n , we always can assume that u n (0) < 
^k~ for all n. Consequently u is not a constant traveling solution of ( 11.3H . □ 



4 ■ 5 Proof of Theorem \2.1\ (v) 



PROOF. Suppose, by contradiction, that for some c G (0,c*), ( II. 3D has a 
traveling wave u(x, t) = u(x + ct) with liminf^_ i . 00 M(£) >> and u(— 00) = 0. 
Thus u(x, t) = u(x + ct) can be larger than a positive vector with arbitrary 
length. It follows from Theorem 12.11 (ii) 

liminf inf u(x,t) > k~ » 0, for < c < c* 

t-KX> \ X \<ct 
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Let c G (c, c*) and x = ct. Then 



lim u{ — (c — c)t) = lim u(—ct,t) > liminf inf u(x,£) >> 0. 

t^oo V ' / t^oo t^oo \x\<tc 



However 



t— >oa 

which is a contradiction. □ 



5 An example 



lim u( — (c — c)t) = u(— oo) = 0, 



Weinberger, Kawasaki and Shigesada [36] established the spreading speed for 
(jl.2p with h{u\) being a unimodal on [0, 1] based on the spreading results for 
cooperative systems in [32]. Our choice of h{u\) is slightly different from [35] 
and simpler than that in [36] . 

Our new contribution to (jl.2p is to characterize the spreading speed as the 
slowest speed of a family of non-constant traveling wave solutions of (11.21) . One 
example of h(u\) in this paper is h{u\) = u\e~ Ul . (jl.2p has two equilibriums 
(0,0), (0, 1) and another coexistence equilibrium. Let U\ = Wi,u 2 = 1 + w 2l 
then (ll.2p can be transformed to 



(5.26) 



dwi r 

-—- = diAwi + w\ Iri - a - 0W\ + riw 2 \ 
at 

dw 2 

— — = d 2 Aw 2 + r 2 (l + w 2 )[-w 2 + h[wi)\ 

In this section, we make the following assumption on h. 

(H4)(i) Assume that h is continuous differentiate on [0, oo) and h(0) = 0, 
h'(0) > 0, h(wi) > 0,tOi 6 (0, oo). Also assume that h m > and h is 
increasing on [0, h m ], decreases on [h m , oo). lim Wl _ >00 h(wi) = 0. 



. i:-,,,,.,.!.. ^,< w -i^ii is strictly decreasing on (0, oo) and lim Wl _ >00 -i^ii 
0. 



(iii) 

h(w 1 ) 2 + 4:h(w 1 )~4:h'(0)w 1 <0, w x e[0,oo). (5.27) 

(H4)(i)(ii) implies that 

h( Wl ) <ti{0) Wl , w x e[0,oo). (5.28) 

We need to verify h{wi) = W\e~ Wl satisfies (H4). h{wi) = Wie~ Wl achieves its 
maximum at h m = 1, and is increasing on [0, h m ] and decreasing on [h m , oo). 
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In addition, h'(0) = 1 and h{w\)/wi = e W1 is decreasing for w\ > 0. It is 
easy to see that e x >x + l,x>0 and e~ x < jtj, a; > 0. Thus, for w% > 

/i(wi) 2 + Ah{w{) - 4ti(0)wi 

< * — + 4wi 

~~ 1W\ + 1 TWi + 1 

iu?(wi + 1) + 4wi(2 Wl + 1) - 4w 1 (2w 1 + l)(wi + 1) (5.29) 

~ (2wi + l)(wi + l) 

-7u>? - 3wj 



(2wi + l)(«;i + i; 



<0 



We also assume that h m < k\ (otherwise, this problem can be dealt as a 
cooperative system) and a < r%. In the nonnegative quadrant, (15.261) has two 
equilibrium (0,0) and (hi, k 2 ) satisfying 

a + Ski = n + nh(ki) 

k 2 = h(h). (5 - 30) 

We claim that (I5.30p has only one positive solution. In fact, the first equation 
of (I5.30p can be rewritten as 

1 = ri ~ a + ri/l(fcl) . (5.31) 

Skt 

From (H4)(ii), ri ~""^ 1 ^ Wl > is strictly decreasing on (0, oo) and 1 = n ~ a ^' 1 ' 1 ' 
has only one solution. 

In order to use Theorem 12.11 we shall define the two monotone systems. As 
indicated in Remark I2.2[ similar ideas for constructing h^ were used in serval 
previous works. However the construction of h~ is different. 



m • 



h(wi), < W\ < h 
h + ( Wl ) [ 

h(hm), Wi > tlm. 

and the corresponding cooperative system is 

— — = diAwi + Wi[ri — a — 8w\ + r 1 w 2 ] 
ot 

— — = d 2 Aw 2 + r 2 (l + w 2 )[-w 2 + h + (wi)] 



(5.32) 



In a similar manner, one can find (I5.32J) has two equilibrium (0, 0) and (kf, k 2 ) 
satisfying 

a + 6kt =n +r 1 h + (kl) 

k+ = h + (k+). (5 ' 33) 
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Fig. 3. The construction of h + and h . The red curve is h. 




Fig. 4. The intersections of n + rih(wi),r\ + r\h (wi) with the line a + 5wi. The 
red curve is r\ + r\h 



Since h + > h, from the first equation of (I5.3ip . it is easily seen that that kf > 
k\. In addition, since k\ > h m , we have k£ = h + (kf) = h(h m ) > h{k\) = k 2 . 

Now there is a h G (0, /i m ] such that /i(/io) = h(ki) and define 



/i (tWi) 



h(wi), < wi < ho, 
h(kf), w\ > ho- 



Then 



< /i - (wi) < /i(wi) < h + ( Wl ) < h'(0)w 1 ,w 1 G (0,^] 
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The corresponding cooperative system for h is 

dwi r 

— — = d 1 Aw 1 + Wi[ri — a — <3u>i + riu> 2 J 

^ 2 (5 ' 34) 

— — = d 2 Aw 2 + r 2 (l + w 2 )[-W2 + h (wi)] 

In a similar manner, one can find (j5.34p has two equilibrium (0, 0) and (fcf , for) 
satisfying 

4=*V). (5 ' 35) 

Similarly, we have fe~ < /ei. In addition, by the definition of /i~, we have 
fc- = /i-(ifcr) < h(kt) < h(h) = k 2 . 



Thus, 

(0,0) « (k^,k 2 ) < (h,k 2 ) < (k+,kt). 

See Remark 12.31 for (Hl)(i). Now it is straightforward to check all other con- 
ditions of (HI) (i)-(iv). 

The spreading results for the cooperative systems were used to establish in 
[36] . We now demonstrate Theorem 12.11 can be used to establish spreading 
speed and traveling wave solutions of the nonmonotone system (I5.26P and 
summarize the results in the following theorem. 

Theorem 5.1 Let d 1 ,a,5,ri,d 2 ,r 2 be all positive numbers. Assume d± > d 2 , 
a < r\, k\ > h m and 

6 > nr2h ' {0) . (5.36) 

r x +r 2 - a 

Then the conclusions of Theorem \2.1\ hold for IJ5.26]) for the minimum speed 
c* = 2y(ri — a)di, A c = ^ c _ 1 Tl > q and v& c , where v\ is defined in 

We now need to check (H2) and (H3). The linearization of (15.261) at the origin 

is 

dwx , . . . 

— — = dxAwi + (n - a)w 1 

d% ( 5 - 3? ) 

d 2 Aw 2 + r 2 (h(0)wi - w 2 ) 



dt 
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The matrix in (jl.8p for (J5.26P is 

d 1 X 2 + r 1 -a 



:? .. ., n \ 



^A = (#) , 

r 2 /i'(0) d 2 A 2 -r 2 ) 



(5.38) 



It is easy to see that (H2) holds. In fact, the principle eigenvalue A\ is ty(A\) = 
d\\ 2 + r\ — a, which is a convex function of A. And 

= gW = rf^ + n-q 
A A 

satisfies the results of Lemma [TTT1 In fact $(A) is also a strictly convex function 



of A. The minimum of $(A) is c* = 2y(n — a)d\. For each A > 0, the positive 
eigenvector of Ax corresponding to \P(A) is 



v\ 



\vl) \ r 2 ft'(0) 



which is also the positive eigenvector of \A\ corresponding to $(A). 
For each c > c*, the left positive solution of $(A) = c is 



c — Jc 2 — Adi (r\ — a) 
c ~ 2d t 

in Lemma 11.11 Further from (15.391) we can see that 

v\ = r 2 h'(0) h'(Q) 

v\ {d\ - d 2 )X 2 + n + r 2 - a a 

where a = 1 + !l=g±&z^ > L 

In order to verify (H3) for (15.321) let 

(w u w 2 ) = (0,0—) » (0,0), e > 0. 
a 

Thus (H3) is equivalent to the following two inequalities 

w\[r\ — a — 8w\ + riw 2 ] < (ri — a)w\ 
r 2 (l + w 2 )[-w 2 + /i + (wi)] < r 2 (ti(0)wi - w 2 ) 

or 



Swi > nw 2 (5.40) 
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and 

ti(Q)wi + w% > h + (w 1 )(l + w 2 ) (5.41) 

It follows that the following equality suffices to verify ( 15.-40]) : 

« > ^^ 

which is true if (15.361) holds. 

In order to verify (I5.4ip . following Weinberger, Kawasaki and Shigesada [36], 
there always is a positive constant <^ such that 

fc'(0) 

w 2 = w 1 . 

2 

Substituting w 2 into (15.411) and multiplying fe ,,Q, at the both sides, we get 

n'[y)wi 

Rearranging the terms produces 

-^ 2 (l - S^") + <^ + K) " h'{0)w! < 0, Wl > 0. (5.42) 

Recall the definition that h + {wi) = h(wi) for w\ < h m and h + {w\) = h(h m ) 
for wi > h m . Since the left side of (15.421) is decreasing in W\ for w\ > h m , we 
only need to verify (I5.42p for W\ <h m . 

By the quadratic formula, the discriminant of (15.421) for W\ < h m is 

h( Wl ) 2 - 4(l - ^L)h f {0) Wl = h( Wl ) 2 + 4h( Wl ) - 4ti(0) Wl , 

which is negative from (I5.27p . Thus the left side of (I5.42p has no real zeros 
and it has to be nonpositive because of the assumption that 1 — h ,i"^' ) > 
for W\ > 0. In fact, (I5.27P is is one of the three possible conditions in [36] to 
guarantee (15.421) holds. 

Notice that h and h are identical around the origin. By the exact same 
arguments (just replacing kf by k~[ and h + by h~), we can verify that (H3) 
holds for ( 15.341) as well. 
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6 Appendix 



In this section, we shall verify provide a direct verification of Lemmas 14.31 
and 14.41 Lower and upper solutions of the equivalent integral equations ( I4.15P 
play a central role in the construction of fixed points of the equivalent integral 
equations through monotone iterations. The lower and upper solutions give 
asymptotic behavior of traveling wave solutions of ( II. 3p . 

Wu and Zou [38|39] . Ma [22~f2"3] verify lower and upper solutions through dif- 
ferential equations, and then use them in monotone iterations of equivalent 
integral equations. While it was pointed out in [3] that the upper and lower 
solutions for differential equations are required to be smooth for delayed equa- 
tions, the author [29] recently directly verified <p + and <p~ are indeed lower and 
upper solutions through the equivalent integral equations for scalar equations, 
where the integrals and compare the two sides of (14.151) were calculated and 
compared. Clearly, in this way, the lower and upper solutions are not required 
to be smooth. 

In this appendix, we shall directly verify that, for rz-dimensional systems, 
<f) + and _ are the lower and upper solutions of (I4.15p . Thus this Appendix 
can be viewed as a continuation of [29] for the direct verification of non- 
smooth upper and lower solutions of the equivalent integral equations for 
n- dimensional systems. 



It should be pointed out that the proof of two lemmas in Ma [221 Lemmas 2.5, 
2.6] can significantly simplify the verification of lower and upper solutions for 
the equivalent integral equations although the conclusions of the two lemmas 
in [22] were about lower and upper solutions for differential equations (see 
[3T| Section 6.1] for more details). As a result, we always can verify them in 
a much simpler way. Nevertheless, a direct verification can provide a further 
evidence that + and 0~ are lower and upper solutions. In addition, by care- 
fully analyzing eigenvalues and corresponding eigenvectors, we identify some 
identities between the parameters and reveal interesting relations between the 
parameters. 

The results in this appendix are natural extensions of those in the author [29] 
for scalar cases. As in [29], for A > let 



(M t (\)) = (3v x -\ 2 Dv x + A x v x 
or 

TV 

Mi(X) = fiu\ - v\d % \ 2 + J2 "K \i = 1, -,N, (6.43) 

i=i 
where v x = {y x ) is the positive eigenvector of jA x in ( 11. 8ft corresponding to 
the principle eigenvalue $(A). For c> c*, recall that $(A C ) = c. It follows that 
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jt c a a c va c = $(A c )z/ Ac = ciy Ac and 

M(A C ) = (J3- dAl + cA c y Ac ,t = 1, ..., JV. 

Because of (JQ21) . Mi (Ac) > 0,i = l,...,iV. Noting M(A) is continuous with 
respect to A, we always can choose a 7 such that 

1< 7 <2, M i ( 7 A c )>0,z = l,...,iV. (6.44) 



In order to simply our proofs, we first prove two identities (Lemmas 16. l\ 
which are the extension of the identities for scalar cases in [21]. Their proofs are 
almost identical to those in [29J except that eigenvector v\ has to be included. 

Lemma 6.1 Assume (HI — H2) hold. Then for each c> c* 

di{Xu + A 2 ») VAij + A c X 2i - A c ) 

PROOF. Recall \ 2l > A u > 2A C > A c . It follows that, i = 1, ..., N, 
M(A C ) f I + 1 



di\A\i + A 2 j) VAij + A c A2i — A 

M(A C ) (Aij + A 2t ) 

di(\\i + A 2 j) AijA2i + (A 2 i — Aij)A c — A c 
Mi (Ac) 1 



di 
M (Ac) 


A _|_ .C. A 

di ^ dt 1Yc 

1 


-A 2 

c 


/3 + cA c - 
<(cA c - 


•rfiA2 
• ^A 2 ) 



(6.46) 



+ cA c - diA\ 



»k c - 



D 

Lemma 6.2 Assume (HI — /72) /io/d and 7 satisfies flff.^j ). T/ien /or eac/i 
c> c*, i = 1,..., A 7 

M(A C ) M(A C ) M( 7 A C ) M( 7 A C ) 

(A w + A> Ac (A 2 i-A c )4 c (Ah + 7 A>; Ac (A 2i - 7 A c )^ Ac ' 



PROOF. Since A 2 j > Ah > 2A C , it follows that A 2 j > 7 A C > A c . Lemma LLU 
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($( 7 A C ) < c) implies that, for i = l,...,N 



M,( 7 A C ) = /?*/ Ac - */ Ac ^( 7 a c ) 2 + E "-' * " 



3 l 3 

7A c a 7 A c 

i=l,...,Af 

US^A.WA.V/. ( 6 - 48 ) 



= (/3-^( 7 AJ 2 + $( 7 A c ) 7 A c )z/ A 
<(/3-rf,( 7 A c ) 2 + c 7 A c )z/; Ac . 

We also note M t (A c ) = ((3 + cA c - <ijA 2 )z/ Ac . Thus, for i = 1, ..., N, we have 
M,(A C )_ + M t (A c ) M 4 ( 7 A C ) M,( 7 A C ) 



(Au + Ac)!/^ (A 2i -A> Ac (A H H- 7 A c )^ Ac (A 2i - 7 A>; Ac 
(Ah + A^M^A,) (\ li + X 2l )M l ( 1 A c ) 



(Ai,A 24 + (A 2l - A H )A C - A2)z/ Ac (AkAm + (A 2l - A l4 ) 7 A c - ( 7 A c ) 2 )z/ Ac 
^^M 4 (A C ) ^±^M,( 7 A C ) 



(| + |A C - A 2 K c (| +| T V-J 7 A c )>> Ac 
v^T4MM 4 (A c ) v^T4MM 4 ( 7 A c ) 



(/3 + cA c - <U 2 K c (/3 + C7 A C - di( 7 A c ) 2 )i/* Ai 
, 2 , ,., / ^(A c ) M i ( 7 A, 



c 



+ P l \ 09 + cA c - d 4 A c 2 R c 09 + c 7 A c - rf,(7A c ) 2 X c ^ 

M,( 7 A C ) 



+ 4/3^(1- 



(/3 + c 7 A c -^( 7 A c ) 2 )z/; Ac 
>0. 

(6.49) 

This completes the proof. □ 



6. 1 Proof of Lemma |^.3| 



The proof of Lemma 14.31 is almost identical to that in [29] for the scalar case 
except that the eigenvector v\ needs to be included, and delay terms are not 
present here. 



In-! 1 



PROOF. Let H = -^, i = 1, ..., N. Then 0+(£) = fe if £ > £*, and 0+(£) = 
v\ e Ac « if % < £*,i = 1,...,N. Note that 0+(£) < v\ e Ac? ,£ e R. In view of 
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(H1)-(H3) we have, for £ < £ 



W + (0) = W(0 + /i(^ + (0) 



1 1 



< M,(A c )e A ^ 



where M»(.) is defined in (I6.43J) . For £ > £*, because of the assumption that 
fll.3p is cooperative and (3 is sufficiently large, we have 



= 0ki. 



Thus, for £ > £*, i = 1, ..., N, we obtain 






oo 



di{X u + X 2i ) J% Jt 



(6.50) 



Thus in view of (14.14[) . we add and subtract the term d . (A ^, A . } jl^, e~ Ali(§ s) <is 
at the left of f l63U|) . Now for i = 1, ...,JV,f > £*, noting that z4 c e Ac ^* = k h 
( 16 .5 Op can be written as 



7-[0 + ](O < h + ,,. \. r (M,(A c ) / ?; e^^e^s 
-Pki f ' e- Xl ^~ s) ds) 






Wj(^li + ^2i) Aij + A c Aij 

fc ie - Alig e A "^ / Mj(A c ) J3_s 

' + di(Aii + A 2i )MAii + A c X ~ W 

* + d,(A u+ S^t)A u ,l ( A "^^ " ^ " ^ 

(6.51) 
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Since Af*(A c ) = (3v{ c - v{ diA 2 c + cA c v\ ,i = 1, ..., N, we have, i = 1, ..., JV 



= <A u (cA c - djA^) - fivj^A 

4/3d, 
: 2^(^ + 4/^ + 
4/3(cA c - d t A 2 c ) - 2(v^+4M + c)/3A 



< ^ 7==^o , ~ (A - <Uc) - M A £ 



v 



A, 



2(y^+4 M + c) (6 . 52) 

. 2(3A c (2c ~ 2d t A c - y/^TI^di - c) 
UAc iy& + 4/3d 4 + c) 



Vk 



2/3A c (c - 2^A C - Vc 2 + 4y9di 



2( Vc 2 + 4/3rf 4 + c) 
<0. 



Combining f !6 . 5 1 1) and (16. 52 p . we see that for £ > £*, z = 1, ..., N, 

7~[0+](O < h. (6.53) 

Similarly, noting z^ c e Ac ^* = ki, one can see that, for £ < £*, i — 1, ..., AT, 
M 4 (A C ) / /■€ 



^ ^/ri t (/ '- m *- ) « a -* 

+ / l e X2 ^- s) e AcS ds) + — — / e X2 ^- s) (3k lS 

Je ' OilAiiH- Ao,: hi 



1 di(\ii + \ 2i ) Je* 

,m i; + a,. 4 



M,(A C ) , e A -g | e A -g 
di(A u + A w ) M M + A c A 2i -A 



e A 2 ^ e -(A 2l -A c )f* x /%. e A2ie e -A 2 i^ 



(6.54) 



A 2 i — A c ' <ij(Aij + \2i) A2i 
^M 4 (A C ) / l l 



+ 



di{\\i + A2i) VAij + A c A2i — A c 



d*(A u + Aai) V(A 2i -A c K c M(A c )A 2i 



Since Mj(A c ) = (3v\ — v\ &iA\ + cA c v\ ,i = l, ..., N, it is easy to see that, by 
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choosing (3 sufficiently large if necessary, 

-k Pkj (-M,(A C ) + ulj)\ 2l - A c uiJ 

; 1 > r i a — T^ — <Gi~ 



(A 2i -A c )z/| c M,(A c )A 2i ' ^{Xv - K)X 2i Mi(A c ) 
= h 



A 2i «A;H - v\K c c) - A c is\f 



u l Kc (\ 2l -A c )\ 2i Mi(K) 
_ h . v i a A 2t (A c ^ -c)-P (6>55) 

1 Ac c <(A 2i -A c )A 2i M,(A c ) 

^ Ac C ^ c (A 2 ,-A c )A 2i M 4 (A c ) 
<0. 

Combining fl6\45|) . fl6\Mj) and fl635|) leads to, for £ < £?, z = 1, ..., JV, 

7^ + ](£) < ^ c e Ac «. 
And therefore, for (gM, 

Wl® <#(0,* = 1.~,W- (6.56) 

This completes the proof of Lemma 14.31 □ 

We now need the following estimate on /, which is an application of the 
Taylor's Theorem for multi-variable functions. Also see 



Lemma 6.3 Assume (H1 — H2) hold. There exist positive constants bij,i,j 
1, ..., N such that 

N N 

fi(u) > J2 d jM°) u j ~J2 b ij( u j) 2 ' u = i u i)iUi e [o,h],i = 1,...,N. 

3=1 3=1 



6. 2 Proof of Lemma \4-4 



Again the proof of Lemma [4.41 is almost identical to that in [2S] for the scalar 
case except that the eigenvector v\ needs to be included and delay terms are 
not present here. 



ln(g^) 



PROOF. Let C = -jtz^T, * = 1, -, N.1f£> £, &"(£) = 0, and for £ < £*, 

4>7(0 = ^ c e A ^-qul rAc e^, l = l,...,N. 
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For £ e R, it follows that 

>0, 
Thus, for £ > £*, 

We now consider the case £ < £*. It is easy to see that 

^ c e Ac « > C(0 > u\ c e^ - g^ Ac e^, £ 6 R,z = 1, ..., JV. (6.57) 

In view of Lemma 16. 3\ ( I6.57p . we have, for £ e R, i = 1, ..., N, 

n n 

> M(A c )e Ac « - gM( 7 A c )e^ - Me 2Ac «, 
where Mj(-) is defined in (16.431) and 

n 

^ = E%KJ 2>0 - (6-59) 

Because of i^ (<£-(£)) > 0, the term /* e A2 ^- s )#;(0-(s))ds of 7I[0"] can be 
ignored in (16.601) . Now we are able to estimate T[0~] for £ < £*,i = 1, ...,N 

WKO > ,,. \ r r ( /* e- A «^- s )M,(A c )e A -d S 

- g / e- AllK - s) M i ( 7 A c )e 7AcS cis 

J — oo 



- Mi / e- A " (? - s) e 2AcS rfs 



+ I * e A ««-^)M i (A c )e AcS ds - g I * e X2 ^- s) Mi(jA c )e^ AcS ds 
- Mi I** e x ^s) e ^s d s 

1 (M t (A c )e A ^ M( 7 A c )e^ 

q- 



di(\ii + \2i)\ Aii + A c ' Aii + 7 A C 

__ e 2A c ^ e A c C*-A 2 ^*+A 2 ^ _ e A c £ 

- Mi- — H M(A C 

Aij + 2A C A c — A2i 

e 7Ac?*-A 2i 5*+A 2 ^ _ e 7A c £ 

- 9 7 r M( 7 A C ) 

7 A C - X 2i 

_ e 2A c C*-A 2 ^*+A 2l ? _ g 2A c ^ 

-Mi 

2A C — A2i 
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(6.60) 



In view of the identity (I6.45p . we subtract two terms to make up a term 
— z/* Ac ge 7Ac ^ and thus we need to add the terms. Recall that 7 A C < 2A C < \ 2 i- 
We ignore two positive terms q-^ .'-,'?,, ,, e (7A e -A 2a )g*+A 2 »£ anc j 



Mi 



(X 2i -2A c )di(Xu+X 2i ) 



e (2A c -A 2i )5*+A 2i 5. Thus. 



'^ HU - d t (X lt + A 2l ) U u + A c + A 24 - A J 6 

Afi(A c ) / 1 1 \z/ 



+ ^ H^^ 



<^(Ah + Am) VAij + A c Am — A c / z/ Ac 
, e^ [ < / M t (A c ) M t (A c ) 

*(A U + Am) | g 7Ac ^(A u + A c K c (Am - A>j 

M 4 ( 7 A C ) Mi( 7 A c 



(6.61) 



(A H + 7A c )^ Ac (A M - 7A c )^ Ac 

g e (2- 7 )A eg _ M t (A c )e(^-^)r c(A2i _ 7Ac)g 
(A u + 2A C ) (A M - A c ) 

M i c (2- 7 )Ac£ 



(A 2i - 2A C ) 



For £ < £*, e^ 2 T)A c ^ e (A 2i 7 A c )g are bounded above. Because of the identity 
f)6.45p . f)6.6ip can be further simplified as, i = 1, ..., N 



m>-](t) > < c ^ - ^w e 



» e l^Z 



e^ A ^ . / M,(A C ) Mi(A c ) 

+ ^(Ah + Am))^^ l(A l4 + A c )z/i c + (A 2i -A>< 



A c 



Mi( 7 A c ) Mi( 7 A c 



(A H + 7 A c )z/ Ac (Am - 7A c )^ Ac j (6.62) 



e (2- 7 )A c £* _ ^M^c^ l e (A 2! - 7 A c )^ 



(Ah + 2A C ) (Am - A 

] ''■ 3 (2-7)A c e* 



(Am - 2A C 



Finally, from (I6.62p and Lemma I6.2[ we conclude that there exists q > 0, 
which is independent of £, such that, for £ < £* and z = 1, ...,N 

W\ (0 > ^A c e Acf - K Ac e 7Ac? . (6.63) 

And therefore, for z = 1, ..., N 

W1(0 >&"(£), eeK- 
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This completes the proof. □ 
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